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Abstract 

We evaluate the Green's function of the D-dimensional relativistic Coulomb 

system via sum over perturbation series which is obtained by expanding the 

exponential containing the potential term V(x) in the path integral into a 

power series. The energy spectra and wave functions are extracted from the 

resulting amplitude. 
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I. INTRODUCTION 



Most physical problems cannot be solved exactly. It is therefore necessary to develop 
approximation procedures which allow us to approach the exact result with appropriate ac- 
curacy. An important approximation method for solving problems in quantum mechanics 
(QM) is the Rayleigh-Schrodinger perturbation theory. It provides us an effective method 
to calculate the approximate solutions of many problems which can not be exactly solved 
by using the Schrodinger equation. Similar to the standard QM, the perturbation method 
can be developed in the path integral framework of QM ]IJ. Historically of utmost impor- 
tance was the application of the perturbation expansion of path integral to the quantum 
electrodynamics by Feynman ||, from which he derived for the first time the "Feynman's 
rules" , which provide an extremely effective method to calculate the perturbation series and 
a clear, neat interpretation of the interaction picture. 

In the past 10 years, perturbation expansion of the path integral has been used to 
obtain the exact Green's functions for (^-function potential problems ]3|-[5],[7|] , non-relativistic 
Coulomb system ||, and to yield the Dirichlet boundary conditions in Refs. ||?] for the non- 
relativistic problems and in Ref. [|K]] for the relativistic problems by summing the 5-function 
perturbation series. 

In this paper, we would like to add a further application of the perturbation method of 
the path integral. We calculate the Green's function of a D-dimensional relativistic Coulomb 
system via summing over the perturbation series. The energy spectra and wave functions 
are extracted from the resulting amplitude. 

II. PATH INTEGRAL FOR THE RELATIVISTIC COULOMB SYSTEM VIA SUM 

OVER THE PERTURBATION SERIES 

Let us first consider a point particle of mass M moving at a relativistic velocity 
in a (D + l)-dimensional Minkowski space with a given electromagnetic field. By using 
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t = —it = —ix^/c, the path integral representation of the Green's function is conveniently 
formulated in a (D + 1)-Euclidean spacetime with the Euclidean metric, 



{.9,J) = diag (l,---,l,c 2 ), 



(2.1) 



and it is given by fTT||T2" 



G(x 6 , x a ; E) = ^ jH dS J Dp® [p] J D D xe- A -l\ 



(2.2) 



The action integral 



dX 



M . 6 . / % ,,. s /A ,(£-V(x)) 2 ,,,Mc 2 

— x< (A) - ! -A(x) • x<(A) - „ (A) 2Mc V' + P(A) — 



(2.3) 



where S" is defined by 
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dXp(X), 



(2.4) 



in which p(A) is an arbitrary dimensionless fluctuating scale variable, and is some 
convenient gauge-fixing functional, such as $ [p] = 5 [p — 1], to fix the value of p(X) to unity 
nyi2]]. h/Mc is the well-known Compton wave length of a particle of mass M, A(x) is the 



vector potential, V(x) is the scalar potential, E is the system energy, and x is the spatial 
part of the (D + 1) vector x = (x, r). This path integral forms the basis for studying 
relativistic potential problems. 

Expanding the potential term V(x) into a power series and interchanging the order of 
integration and summation, we obtain the result 

%n r°° .„ r , -*/£dAp ( A)£ 



G(x 6 , x a ; E) = ^ fc | o °° J Dp® [p] eT* A. 1 dAp(A) ^(x fe , x a ; A b - A a ) 



2Mc Jo 

with the series expansion of the pseudotime propagator 



(2.5) 



K(x 6 ,x a ; X b - X a ) 



n=l 



x / D u xe 



^x' 2 (A)-itA(x).x'(A)-p(A)fI^ 
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x / Ab dA lP (A 1 )V r (x(A 1 )) [ Xb dA 2 p(A 2 )y(x(A 2 )) • • • dA n p(A n )V(x(A n )) 



(2.6) 



where we have defined the quantities (3 = E/Mc , £ = ( M 2 c — E 2 )/2Mc , and 



Ko(xfe,x a ; A fe - A a ) 



n -rf b d\ 
D D xe Rha 



^L T x^(A)- J fA(x).x'(A)-p(A)|^ 



(2.7) 



Ordering the A as Ai < A 2 < • • ■ < A n < A and denoting x(Afc) = X&, the perturbative series 
in Eq. ( |2.6| ) turns into M 



K(x b , x a ; A fe - A a ) = K {x b , x a ; X b - A a ) + > — r / ^ / dA n _i ■ • ■ / dAi 

^ \ ft / JA a JA a JA a 



X 



II ^o(xj+i, Xj-, Xj+i - Xj) 

j=0 



Y[ PkV(x k )dx k , 



(2.8) 



k=l 



where A = A a , A n+ i 
potential, we have 



A ,x„ + i = x&, and x = x a . In the case of an attractive Coulomb 



A(x) =0, V{r) 



(2.9) 



The perturbative expansion in Eq. (|2.8|) becomes 

~ ( (3e 2 \ n f x b /-An fAa 

K(x 6 , x a ; A - A a ) = ^ (x o , x a ; A - A a ) + > — / d\ n / c/A n _i • • • / d\i 

~i \ Ti J Jx a Jx a J\ a 



H -^o( x i+i) x jn ^i+i — Xj 

j=0 



The corresponding amplitude K takes the form 



k=l 



r k 



K (x o ,x a ; A 6 - A„) 



D D xe hJxa 



2p(A) 



x' (A)-p(A), 



(2.10) 



(2.11) 



where a = e 2 /he is the fine structure constant. We now choose $ [p] 
value of p(A) to unity. The Green's function in Eq. (|2.5|) becomes 



5 [p — 1] to fix the 



G(x& ' x - E) = 2k/ 



dSe i 



K (x h ,x a ; S) 



^ / Be 2 \ n r x b /-a„ r x 2 r 

~l \ ft J JXa JXa JXa J 



i=o 



n — 

fc=l rfc 



> . 



(2.12) 



We observe that the integration over S is a Laplace transformation. Because of the convo- 
lution property of the Laplace transformation, we obtain 

ih 

G(xfe,x a ;E) 



x 



2Mc 



fc=i ' k 



(2.13) 



We now perform the angular decomposition of Eq. fl2.13| ) [|12| — 13~4 1 . This can be reached by 
inserting in Eq. (|2.13| ) the expansion of Go i n term of the D-dimensional hyper spherical 
harmonics Yj m (x) ||15|| : 

M 



Co(x J+ i, Xj-; S) 



TE^(^ + i^,^)E^m(x J+1 )^(x,), (2.14) 



where the gf is given by |T4 



^ -fS -M(rj +1 +rf)/2hS t 

o S 6 6 V('+^-i) a - a V ft 5 



fMr j+irj 



(2.15) 



The notation / denotes the modified Bessel function. Integrating over the intermediate 
angular part of Eq. (|2.13|) , we arrive at 

ih 



G(x b ,x a ;£) = __^G z (r b ,r a ^)5:il m (x b )F4(x a ). 



(2.16) 



z=o 



The pure radial amplitude Gi(r b ,r a ; £) has the form 

M 1 ^ (M(3e 2 ^ n 



Gi(r b} r a ;S) 



h (r b r a ) D /*- 



n=0 



9i n \ r b,r a ;S) 



h 2 



(2.17) 



with <7j is given by 



9i l '(r b ,r a ;£) 



oo roo 



J 



n^ (0) ( r i+i' r i^) 



n rfr *=- 
/c=i 



f2.18) 
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To obtain the explicit result of gj n \ we note that 

_ e .fS e -M(r 6 2 +r2)/27i.Sj 



00 dS s , 



S 



2 / cfe— _ e - K ^+^) cothz /„ 



(Mr b r a 

y/(l+D/2-l)*-a* \ % s 



(2.19) 



/o ~~sinh2~ 2 v /(/+^/2-i) 2 - Q 2 ^ ginhz ) 

with k = VMV - E 2 /hc. The equality in Eq. (gl9|) can be easily proved by the formulas 



e 2i ^ 

— ; — exp 
smh y 



--((a + (b)cothy 



sinhy 



r((i + /x) /2-1/) r , , , 
^r(, + i) ^/2(*C 6 )M^ /2 (tc a ) 



with the range of validity 

Cb > (a > 0, 

Re[(l + n)/2 -u]>0, 
Re{t) > 0, 1 argt |< ir, 

where M^ u and W^ v are the Whittaker functions, and 

,2ab 

y 



r %-» e -(^)/v Iv ^) = 2I u (2aVz)K v {2by/z), 
Jo y y 

with the range of validity 



a < b, Rez > 0. 



From Eq. ( [2.19 ), we get, by using the formula 



the result 



& (1) (r b ,r ;£) = I gf > {r b) r- ) £)gf , {r,r a - ) £)dr 
Jo 



M, 



2^ PCX) 

— / zh(z)dz, 

K Jo 



^2.20) 



(2.21) 



[2.22) 



(2.23) 



where the function h(z) is defined as 

1 



/,, - ) - ■ -K(r h +r a )cathz T ( 2K^r b ^\ l 

{ ) shlhz ^^/(l+D^-l) 2 -* 2 I sinh ^ 



(2.24) 



The expression for g\ n \r b , r a ; 8) can be obtained by induction with respect to n, and is given 
by 



g\ n \r bl r a -£) 



2 n + 1 X roo 



n\ k Jo 

Inserting the expression in Eq. ( |2.17| ), we obtain 

Gi(r b ,r a ;£) 



z n h(z)dz. 



(2.25) 



ft {r b r a ) D/2 ~ 



poo 

x / dze 
Jo 



( 2Mf3e 2 \ . 



1 -K(r b +r a )cothzT ( 2K^/r b r a 

sinhz 2^(2+^/2-1)2-^ I sinhz 



[2.26) 



With help of the formula in Eq. ( [2.20| ) , we complete the integration of Eq. ( p.26| ) , and find 
the radial Green's function for r b > r a in the closed form, 

1 Mc 



Gi(r b ,r a ;E) 



( n r a )(D-i)/2 ^/ M 2 c a _ E 2 



x ■ 



r 1/2 + J(/ + d/2 - 1) 2 - « 2 



Ea 



VM 2 c*-E 2 



r (i + 2^/(1 + d/2 - r 



2 - a 2 



xW 



\ / M 2 c 4 -E 



^^ { l+D/2-l) 2 - a 2 { hc 



—VM 2 c A - E 2 



n 



(2.27) 



The energy spectra and wave functions can be extracted from the poles of Eq. Q2.27I) . 
For convenience, we define the following variables 



k = ±VM 2 c 4 - E 2 , 



v 
I 



he 

aE 

VM 2 c 4 -E 2 ' 



(2.28) 



(Z + D/2 - l) 2 - a 2 - 1/2, 



From the poles of Gi(rb, r a ; E), we find that the energy levels must satisfy the equality 



— v + I + 1 = — n r , n r = 0, 1, 2, 3, 



(2.29) 



Expanding this equation into powers of a, we get 



E nl « ±Mc 2 < 1 - - 



1 2 



71 + 1/2(^-3) 



a 



[n + l/2(D-3)] ; 



x 



1 



3 



1 



+ 0(a e 



^2.30) 



_2 [I + 1/2(D - 2)] 8 [n + 1/2(D - 3)] _ 
Here n is defined by n r — n — I — 1. We point out that by setting D=3, the energy levels 
reduce to the well-known form 



+ 0(a 6 ) . 



C2-31) 



.2Z + 1 8nJ 

The pole positions, which satisfy u = fii = n+l—l (n = l+l, 1+2, 1+3, ■ ■ •), correspond to 
the bound states of the D-dimensional relativistic Coulomb system. Near the positive-energy 
poles, we use the behavior for v n\, 



~ M (-) n r 1 / E \ 2 2hMc 2 



£ 2 - K 



(2.32) 



with dn = a-R^jf- being the modified energy-dependent Bohr radius and n r = n — I — 1 
the radial quantum number, to extract the wave functions of the D-dimensional Coulomb 

system 



Gi(r b ,r a ;E) 



( rb r a ){D-i)/2 n 



ir +1 \Mc 2 J 



E \ 2 2hMc 2 



E 2 - El 



x- 



i 1 (ni + l)\ __ (n+rt 



(2l + l)\ 



l2 ^2 



nfa H (n — I — 1)! 



H-i 



2ty ~ 2r 

xM(-n + 1 + 1,21 + 2; -)M(-n + 1 + 1,21 + 2; - 

a H ni a H ni ' 



i ~ , E sQ 2/iMc 2 



(— — ) 



(2.33) 
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where we have expressed the Whittaker function M\ tfJ ,(z) in terms of the Kummer functions 
M(a,b; z), 



M x ,„(z) = ^ +1/2 e" 2/2 M(/i - A + 1/2, 2/i + 1; z). 



(2.34) 



From this we obtain the radial wave functions 
Rniir 



1 



(ni + l)\ 

nfa 1 ^ (2[+l)!\|(n-Z-l)! 



i+i 



x ( —^-) e- r/ ' aHhl 
\a H niJ 



2r 

M(-n + l + l,2l + 2; ' 

a H ni' 



(2.35) 



The normalized wave functions are given by 



1 



^nzm(x) = {D _ l)/2 R n i(r)Y lm (±). 



(2.36) 



Before extracting the continuous wave function we note that the parameter k is real for 
E |< Mc 2 . For | E |> Mc 2 , the square root in Eq. (|2.28|) has two imaginary solutions 



1 



k = Tilt, k = -VF^MV, 
Tic 



(2.37) 



corresponding to 



v = v 



Ea 
hck 



[2.38) 



Therefore the amplitude has a right-handed cut for E > Mc 2 and E < —Mc 2 . For simplicity, 
we will only consider the positive energy cut. 

The continuous wave function are recovered from the discontinuity of the amplitudes 
Gi(rb,r a ; E) across the cut in the complex E plane. Hence we have 

discG;(r b , r a \ E > Mc 2 ) = G t (r b , r a \ E + irj) - Gi(r b , r a \ E — irf) - 



( r6 r a )P-D/2 



M 



-iv 



l + l) 



(2l + l)\ 



^r + i/ 2 (-2^r 6 )M i5j[+1/2 (-2^r a ) + 



v — > —v) 



^2.39) 



Using the relations 



M K ,,(z) = e^^M^i-z), (2.40) 
where the sign is positive or negative depending on whether Imz > or Imz < 0, and 



Wx Jz) = e ftrA e -«r(H-l/2) r ^ + A+l/2) 



X 



r(2/i + 1) 



(2.41) 



r(/i- a + 1/2) 

which is valid only for arg (z) E (—ir/2, 3ir/2) and 2/i ^ —1, —2, —3, ■ • ■. The discontinuity 
of the amplitude is found to be 



discdin, r a ;E > Mc 2 ) 



i M 


r 


'-iu + l + l) 


2 


{r b r a ){D-i)/2 h ~ k 


r(2/~+2) 


2 



xe™M_ lDj+l/2 (2tkr b )M lDj+l/2 (-2tkr a ). 



(2.42) 



Thus we have 



/ 

J A 



discGi(r b ,r a ;E > Mc 2 ) 

discG/(r b ,r a ;£ > Mc 2 ) 



Mc 2 27r/l 

1 /■«> (ftc) 2 M; 



2vr^ 7-oo ^m^ 4 + (ftcA;) 2 

2 /-oo _ ^ 

^(^)^(n)^(r a ). 



(2.43) 



From this, we obtain the continuous radial wave function of the D-dimensional relativistic 
Coulomb system 

T~ i |r(-«/ + z + i)' 




1 , ( chk \2 



1/2 




(2Z + 1)! 

r(-^ + z~+ 1) 



e^ /2 M^~ +1/2 (-2^r) (2.44) 



1 , / cfifc 
1 ^ 



1/2 



(2Z + 1)! 

x e wi>/2 e ifcr (-2iAx) ,+1 x M(-ii/ + Z~+ 1, 2/ + 2; -2ifcr). 



(2.45) 



It is easy to check the result is in accordance with the non-relativistic wave function when 
we take the non-relativistic limit. 



10 



III. CONCLUDING REMARKS 



In this paper we have calculated the Green's function of the relativistic Coulomb 
system via sum over perturbation series. From the resulting amplitude, the energy levels 
and wave functions are given. Different from the conventional treatment in path integral 



using the space-time and Kustaanheimo-Stiefel transformation techniques (e.g. 



the 



method presented here just involves the computation of the expectation value of the moments 
Q n (Q = dAp(A)V(x(A))) over the measure 



and summing them in accordance with the Feynman-Kac formula [17 



G(x fc , x a ; E) = —j o dSj Dp® [p] e 



x E 



exp 



l r Xb 



d\ P (X)V{x(X)) 



%h 



2Mc Jo 



poo p 

j o dS J Dp® [p] e 



(3.1) 



x±tm n E 



n=l 



nl 



d\p(X)V(x(X)) 



(3.2) 



where the notation E [★] stands for the expectation value of the moment *. 

We hope that the procedure presented in this article may help us to obtain the results 
of other interesting relativistic systems. 
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